First-principles investigation of effect of pressure on BaFe 2 As 2 
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On experimental side, BaFe2As2 without doping has been made superconducting by applying 
appropriate pressure (2-6 GPa). Here, we use a full-potential linear augmented plane wave method 
within the density-functional theory to investigate the effect of pressure on its crystal structure, 
magnetic order, and electronic structure. Our calculations show that the striped antiferromagnetic 
order observed in experiment is stable against pressure up to 13 GPa. Calculated antiferromagnetic 
lattice parameters are in good agreements with experimental data, while calculations with nonmag- 
netic state underestimate Fe-As bond length and c-axis lattice constant. The effects of pressure 
on crystal structure and electronic structure are investigated for both the antiferromagnetic state 
and the nonmagnetic one. We find that the compressibility of the antiferromagnetic state is quite 
isotropic up to about 6.4 GPa. With increasing pressure, the FeAs4 tetrahedra is hardly distorted. 
We observe a transition of Fermi surface topology in the striped antiferromagnetic state when the 
compression of volume is beyond 8% (or pressure 6 GPa), which corresponds to a large change of 
c/a ratio. These first-principles results should be useful to understanding the antiferromagnetism 
and electronic states in the FeAs-based materials, and may have some useful implications to the 
superconductivity. 

PACS numbers: 



I. INTRODUCTION 

The discovery of high temperature superconductivity 
in iron-pnictide materials^ has attracted widespread in- 
terest in elucidating the superconductivity mechanism 
and searching for new high-T c materials^^. The high- 
est critical temperature (T c ) observed in the RFeAsO se- 
ries compounds (R= La-Gd) is up to 55 K with electron 
doping^. Recently, superconductivity was discovered in 
the AFC2AS2 materials (A = Ca, Ba, Sr) which have sim- 
ilar FeAs layers. It has been reported that the transition 
temperature can be as high as 38 K for (K,Ba)Fe2 AS2 and 
(K,Sr)Fe2As2^. All of these compounds show an antifer- 
romagnetic order and a structural transition with similar 
ordering temperatures^. Although the superconductivity 
mechanism is not yet known up to now, the phonon me- 
diated pairing seems to be ruled outi£. The fact that an- 
tiferromagnetic spin fluctuations coexist with supercon- 
ductivity indicates that the nature of superconductivity 
should be complex and unconventional in the FeAs-based 
compounds&^iii, like that of the cuprate high-T c super- 
conductors. 

Stoichiometric BaFe2As2 is not superconductor but an 
antiferromagne t 12 ' 13 . Superconductivity can be achieved 
by cither electron or hole doping. In hole doped super- 
conducting (K,Ba)Fc2As2, T C =38K has been obtained by 
replacing 40% Ba by K£. So that BaFe 2 As 2 is a typical 
parent compound of the high temperature iron-arsenides 
superconductors. In addition, recent experiments con- 
firmed that BaFe2As2 without doping becomes supercon- 
ducting under high pressure^. The pressure inducing 



superconductivity was also observed in the same group 
compounds SrFe2As2^ and CaFe2As2^. It is very inter- 
esting to investigate the properties of AFe2As2 materials 
under high pressure. 

In this paper we use full-potential augmented-plane- 
wave method based on the density-functional theory to 
investigate pressure effect on the crystal structure, mag- 
netism and electronic structure of BaFe2As2- The rea- 
son we concentrate on BaFc2As2 to investigate the pres- 
sure effect is that its superconductivity, in a wide pres- 
sure range (between 2-6 GPa) and the highest T c up 
to 28 K, is more robust than that of the Sr and Ca 
compounds 1 —. At first, we analyze the effects of differ- 
ent magnetic orders on the equilibrium lattice parameters 
and corresponding electronic structures. Our calculation 
indicates that the striped antiferromagnetic state is the 
ground state, its calculated crystal structure and lattice 
constants are consistent with experimental data. With 
increasing pressure, there is a topological transition of 
Fermi surfaces in the striped antiferromagnetic state ap- 
proximately at 6 GPa. This Fermi surface transition is 
possibly related to the structural transition and consis- 
tent with the loss of superconductivity at the pressure. 
These full-potential DFT results should be helpful to un- 
derstand the pressure-induced properties of BaFe2As2^i. 

The remaining part of this paper is organized as fol- 
lows. In next section, we describe our computational 
method and parameters. In section III, we analyze the 
crystal structure, magnetic orders and electronic struc- 
tures under ambient pressure. In Section IV, we present 
pressure-driven changes of the crystal structure, mag- 
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netic order, and electronic structures. In section V, we 
compare our calculated results with known experimental 
results. Finally, we give our conclusion in section VI. 



II. COMPUTATIONAL METHOD 

The present calculations are performed by using an 
ab initio all-electron full-potential linearized augmented 
plane wave (FLAPW) method based on the density- 
functional theory, as implemented in WIEN2K codei^. 
Both local density approximation (LDA)i£ and general- 
ized gradient approximation (GGA)i£ are used in our 
calculations. The atomic sphere radii of 2.4, 2.0, and 
1.9 atomic unit are used for Ba, Fe and As, respec- 
tively. The self-consistent calculations are considered to 
be converged only when the integrated charge difference 
between input and output charge density is less than 
0.0001. The internal parameter za s , describing the As 
position, is relaxed until the force per atom is smaller 
than lmRy/a.u. We use 500 /c-points for Brillouin-zone 
integrations in the antiferromagnetic calculations, and at 
least 1500 Appoints in the nonmagnetic calculation. 

BaFc2As2 has ThCr2Si2 type body centered tetrago- 
nal structure (/4/mram) with lattice parameters a = 
3.9625 A and b = 13.0168 A at room temperature. The 
FeAs layers are stacked so that the As atoms face each 
other and the Ba atoms sit in the resulting 8-fold coor- 
dinated square prismatic sites between them. At 140 K, 
it undergoes a structural phase transition from tetrag- 
onal to orthorhombic structure (Framm) with lattice 
constants a = 5.6146 A, b = 5.5742 A, and c=12.9453 
A (at 20K), accompanied by a magnetic phase transi- 
tion from a Pauli paramagnetism to an antiferromag- 
netic order. Neutron scattering experiments have shown 
that the magnetic moments are aligned ferromagnetically 
along 6-axis but antiferromagnetically along a-axis and c- 
axi s 13 i 19 , the named striped antiferromagnetic order state 
(AF1). We also calculated checkerboard antiferromag- 
netic order state (AF2), in which magnetic moment of 
nearest-neighbor Fe are antifcrromagnctical order while 
magnetic moment of next nearest-neighbor are ferromag- 
netical order. In all the calculations, we keep b/a ratio 
in term of the experimental value because we have found 
that the changing of b/a has tiny influence on the total 
energy and electronic structure. 



III. ANALYSIS OF AMBIENT-PRESSURE 
PROPERTIES 

As a starting point, we took the experimental lat- 
tice parameters of tetragonal paramagnetic state and the 
striped antiferromagnetic state of BaFe2 As2 , and relaxed 
the internal coordinate of As by LDA calculation with 
generalized gradient approximation (GGA). For nonmag- 
netic state, we found that relaxed internal parameter of 
ZAs=0.3448 is much lower than the experimental value of 
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FIG. 1: Band structure (top two panels) and Fermi surface 
sections (the other four panels) of BaFe2As2 of experimental 
lattice constant with experimental As coordinates za s (left 
column, labelled as exp. str) and with ab initio optimized 
value (right column, labelled as opt. str). The middle pan- 
els show the Fermi surface sections extended to the second 
Brillouin zone in the k x direction for the k z =0 plane, and the 
bottom panels show those for the k z =n /c plane. 



ZAs=0.3545. While for the striped antiferromagnetic or- 
der, we arranged the direction of iron magnetic moment 
in terms of experimental observation^, and got a value 
of za s =0.3515, which is in better agreement with exper- 
imental data of zas=0. 3538^2. Other antiferromagnetic 
settings with optimized lattice parameters cause smaller 
ZAs value, but are better than nonmagnetic result. Calcu- 
lated magnetic moment of iron inside muffin tin sphere 
is 1.86/xb, almost twice larger than the value of 0.8/is 
measured by neutron scattering experiment^. 

It turns out that the band character of BaFc2As2 is 
similar to that of LaFeAsOiS, except that La / and O 
p staes are missing in BaFc2As2 compound. The band 
splitting of ds z 2_ 1 states is also significantly weaker in 
BaFe2As2 than in LaFcAsO, which indicates that inter- 
layer bonding through the Ba layer of BaFe2 As2 is weaker 
than the LaO layers of LaFcAsO. Since there have been 
many reports on the electronic structure a 21 ' 22 i 23 ' 24 , we 
show only the bands near Fermi level to manifest the key 
parts. The band structure and Fermi surfaces of tetrago- 
nal phase are illustrated in Fig. [1] for both experimental 
and optimized structures. 

Fermi surfaces of both experimental and optimized 
structures have a doubly-degenerated electron pocket 
centered at the M point. These sheets have a domi- 
nant d xz ,d yz character, only a little component of d xy 
character. Around T point, there is significant differ- 
ence. For optimized structure, Fermi surfaces comprise 
a doubly-degenerated cylindrical hole pocket centered at 
the r point, mainly consisting of d xz i yz character and lit- 
tle mixture of d x 2_ y 2 and d 3z 2_ 1 . While for experimen- 
tal structure, another new cylindrical hole sheet, mainly 
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FIG. 2: (color online). Band structures and Fermi surfaces of 
the striped antiferromagnetic state with experimental lattice. 
Here, we use a orthogonal lattice and show Fermi surfaces 
on a conventional orthogonal BZ. The choice of BZ does not 
influence our main discussions and results, but is better to 
guide for eyes. 



d xy character, inserted to the two hole sheets. Those 
sheets have more two-dimensional character. Compared 
to the observation of APRES experiment^, the calcu- 
lated results using the experimental lattice parameters 
are in better agreement with ARPES data. Furthermore, 
in comparison with other calculated result a 21 ' 22 i 23 ' 24 , we 
would find that whether d xy bands cross Fermi level or 
not, it is sensitive to za s , which is similar to the results 
discussed in LaFeAsO syste m 11 ! 23 . 

Furthermore, it is found that the band width of the d xy 
bands changes little, but they shift downwards by about 
0.2 eV in comparison to the d xz / yz bands. This fact indi- 
cates that the tiny change of As position has a significant 
effect on the crystal splitting between d xy and d xz / yz or- 
bital bands. Since there is direct d xy coupling between 
the two Fe atoms, the d xy bands are not sensitive to the 
change of the za s position and are more likely affected 
by the Fe-Fe distance, while d xz / yz orbital are likely cou- 
pled with As p orbital states (mainly p x and p y ), the 
dxz/yz bands have stronger k z dispersion than the d xy 
bands and the corresponding Fermi surfaces have more 
three-dimensional (3-D) character, which change more 
distinctively as the Fe-As bond length is shortened. The 
discrepancy of Fe-As bond length between ab initio cal- 
culations and experimental data, which is induced by the 
magnetism, is a common feature of superconducting iron- 
pnictide family. This is because the iron-pnictide mate- 
rials is in proximity to a quantum critical point, results 
in a strong spin fluctuations^. 

In Fig. [21 we show the band structures and Fermi 
surfaces of the striped antiferromagnetic state with ex- 
perimental lattice constants. In Fig. [31 we show the 
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FIG. 3: (color online). Top panel: total DOS for the striped 
antiferromagnetic state. Bottom panel: spin resolved Fe 3d 
DOS, showing majority (blue) and minority (red), and the As 
4p DOS (green). 



total density of states (DOS) and the partial DOS of Fe 
3d and As Ap states in different boxes. The bands from 
~ — 6 cV to —2 cV arc mainly from As p states. The Fe 
d bands range from —2 cV to +3 cV, being above As p 
states. Around —3 cV, there is significant hybridization 
between the Fe d states and As p states. The band char- 
acter around the Fermi level is dominated by d orbital. 
Due to the striped antiferromagnetic order, a pseudo- 
gap is opened near the Fermi level around T point. In 
Fig. [21 we could find the shape of Fermi surfaces is not 
regular due to the complex d bands dispersion around 
the Fermi level induced by antiferromagnetic interaction. 
There are hole-type Fermi surfaces around Z point (yel- 
low), two slender electron- type Fermi surfaces parallel 
to T-Z line located at k x direction, and two irregular 
electron-type Fermi surfaces located at T-Yl line (pink) 
which arc separated and confined in one quarter of BZ. 
They are mainly of d xy , d xz and d yz characters, with 
some mixing of d x 2_ y 2 and c? 3z 2_ 1 characters. 

Compared to the quasi-two-dimensional Fermi surfaces 
of the striped antiferromagnetic LaOFeAs^, the Fermi 
surfaces of BaFe2As2 have distinct three-dimensional 
topology. The effect of the orthorhombic lattice dis- 
tortion on the electronic band structure as well as the 
Fe moments is very weak. The intcr-layer spin interac- 
tion is antiferromagnetic because our total energy calcu- 
lations show that the ferromagnetic order state is higher 
by about 2.6 meV in total energy per Fe atom. The inter- 
layer spin coupling constant is estimated to be J± = 0.8 
meV. The most stable striped antiferromagnetic order 
state with a magnetic moment is energy lower 124 

meV per Fe atom than the nonmagnetic state, while the 
checkerboard antiferromagnetic order state with a mag- 
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FIG. 4: Calculated values of total energy (a), pressure (b), 
c/a ratio (c), As-Fe-As angle (d), zas (e), Fe-As bond length 
(f), As-As bond length (g), Ba-As bond length (h), Fe-Fe 
bond length (i), magnetic moment of iron (j) as functions of 
volume (Vb is the experimental volume of the antiferromag- 
netic state per formula unit, which is about 0.2% smaller than 
the experimental value of nonmagnetic state). Filled symbols 
indicate the striped antiferromagnetic results, and open sym- 
bols nonmagnetic ones. For the As-As distance and As-Fe-As 
angle, due to the orthorhombic distortion, antiferromagnetic 
state has three values while nonmagnetic state has only two. 



netic moment 1.63/kb is only 54 meV energy lower. By 
mapping to a simple Heisenberg model only including the 
1st and 2nd neighbor exchange couplings, we could esti- 
mate the coupling constants: J\ a = 28.2 meV, Jit = 27.3 
meV and J2 = 17.4 meV. The tiny difference between the 
coupling constants J\ a (along a-axis) and Ju, (along b- 
axis) indicates that exchange coupling of BaFe2As2 is not 
so anisotropic. 



IV. EFFECT OF PRESSURE ON CRYSTAL 
STRUCTURE, MAGNETIC ORDER, AND 
ELECTRONIC STRUCTURE 

Now we investigate the effect of pressure on the crys- 
tal structure and electronic structure of BaFe2As2- We 
compress crystal volume up to 18%, corresponding to 
pressure of 18 GPa (about 180 kbar). We relax the in- 
ternal parameters zas and c/a ratio (keeping b/a ratio 
of orthorhombic structure in term of the experimental 



TABLE I: Calculated equilibrium lattice constant (in unit of 
A)with internal parameter (zab) of BaFe2As2 obtained from 
the fully optimized nonmagnetic (NM) and striped antiferro- 
magnetic (AF1) calculations. 

exp.(NM) exp.(AFl) cal.(NM) cal.(AFl) 



Spacegroup 

o(A) 

6(A) 

C (A) 

Z As 

rf(As-Fe)(A) 



M/mmm 
3.9625 

13.0168 
0.3545 
2.403 



Fmmm 
5.6146 
5.5742 

12.9453 
0.3538 
2.392 



I4/mmm 
3.9591 

12.6210 
0.3448 
2.3330 



Fmmm 
5.6378 
5.5931 
12.9624 
0.3512 
2.3789 



value) of all structures at intervals of 2% volume com- 
pression. We subsequently fit the E~V curves to a Birch- 
Murnagham equation-of-state, to obtain the equilibrium 
volume Vo and the bulk modulus Bq for each system. 
The magnetic moment decreases with increasing pres- 
sure, which results in possible small change of b/a, but 
this change hardly influences other structural parameters 
such as zas and the electronic structure. Therefore, it is 
reasonable to keep b/a in our magnetic calculations. 

The theoretical lattice constants and internal parame- 
ters obtained by the complete optimization are listed in 
Table [U As expected, nonmagnetic calculations underes- 
timate c/a ratio because As atom is more close to Fe lay- 
ers. We could find that nonmagnetic calculation roughly 
predict the in-plane Fc-Fc bonding correctly, but shorten 
out-of-plane Fe-As bond length, therefore As atom is 
more close to iron layer, resulting in a small zas- While 
the calculated lattice parameters of the striped antiferro- 
magnetic state are in better agreement with experimental 
data. Furthermore, the striped antiferromagnetic calcu- 
lation could also describe out-of-planc Ba-As and As-As 
bond better. 

In Fig. |4l we show the dependence of structural param- 
eters on the compression of volume for both nonmagnetic 
and striped antiferromagnetic states. The total energy, 
pressure, the bond length, bond angle, and the magnetic 
moment of iron are illustrated as functions of the com- 
pression of volume. The calculations indicate that the 
striped antiferromagnetic state is the ground state as long 
as the pressure is less than 15 GPa, about 16% compres- 
sion of volume, then the spin-polarization disappears. It 
is found that ferromagnetic state is unstable under a large 
range of pressure. The checkerboard antiferromagnetic 
order is always higher in total energy due to the smaller 
calculated magnetic moment, which results in smaller 
magnetic exchange energy It also could not give reason- 
able lattice parameters, for example, to predict smaller 
c/a ratio and za s - As the volume is compressed by 8%, 
the corresponding magnetic moment of the checkerboard 
antiferromagnetic order state decrease from 1.6 [is to 
0.5 hb- After having tried several other antiferromag- 
netic configurations, we find that the striped antiferro- 
magnetic state is the most stable due to its largest mag- 
netic moment. Some of other configurations even could 
not converge to magnetic solutions. 
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The high pressure experiment indicates that supercon- 
ductivity appears at about 2.5 GPa, then T c reaches a 
maximum near 3.5 GPa, and disappears above 5.6 GPaM. 
As shown in Fig. 0Jb), the compression of volume could 
be estimated to cover such a range of pressure. In terms 
of the antiferromagnetic results, as volume is compressed 
from 2% to 8%, the pressure range relevant to supercon- 
ducting dome is covered completely, while in the non- 
magnetic calculation, the volume needs to be compressed 
more, from 4% to 10%. 

It is clear that Fe-As bond length is enhanced signifi- 
cantly than Fe-Fe after considering the spin-polarization. 
Hence, As atom goes away from Fe layer giving a larger 
zas- With a 10% compression of volume, the magnetic 
moment of iron deceases from 1.8 [1b to 1.2 /is, which 
reduces the discrepancy of Fe-As bond length between 
nonmagnetic and antiferromagnetic states. In contrast 
to tendency of iron related bonding, Ba-As bond length 
becomes shorter after considering the spin-polarization. 
In the striped antiferromagnetic state, chemical bonding 
becomes more rigid against pressure, therefore c/a ra- 
tio changes less than nonmagnetic state, as a result, the 
compressibility is less anisotropic. 

Although the Fe-Fe bond length and Fe-As bond length 
decrease significantly with the increase of pressure, the 
As-Fe-As bond angles e change slightly. Under 10% 
compression of volume, three As-Fe-As angles of anti- 
ferromagnetic lattice are found to be 107.6°, 108.3° and 
112.5°, roughly changed one degree compared to the am- 
bient condition. Former two values correspond to the 
As-Fe-As angles whose two As atoms separated by the 
iron plane, and the latter one represents two As atoms 
on the same side. Due to the stronger Fe-As bonding 
in the antiferromagnetic state, the FeAs4 tetrahedra gets 
more elongated along c-axis, so that FeAs4 tetrahedra in 
the antiferromagnetic state is distorted less than that in 
nonmagnetic state, but it is still far away from the ideal 
structure. 

In Fig. [S]we show the band structures and Fermi sur- 
faces with 2%, 4%, 6%, 8% and 10% compression of vol- 
ume. We can find that the Fermi surfaces expand signif- 
icantly due to the expansion of band width by pressure. 
As volume compressed to 4%, corresponding to 3.1 GPa 
pressure where superconductivity has been observed, the 
Fermi surfaces lie on antinodal directions of V - Yl line 
expand and across T - M line to invade neighbor quarters 
of BZ. Thus, when pressure increases to 2.5 GPa where 
the superconductivity takes place, the main change of 
Fermi surfaces of the striped antiferromagnetic state is to 
expand in BZ. Furthermore, as volume compressed over 
8%, the Fermi surfaces centered around Z point across T 
point to form a connected cylindrical tube. The topolog- 
ical change of Fermi surfaces should respond to the slope 
change of c/a ratio about 8% compression of volume (see 
Fig. [4j. It is notable that the superconductivity disap- 
pears as pressure is higher than 6 GPa, that is roughly 
same with the place where the Fermi surfaces topological 
transition occurs. 




FIG. 5: (color online). Band structures (left column) and 
Fermi surfaces (right column) of the striped antiferromagnetic 
state with 2%, 4%, 6%, 8% and 10% (from top to bottom) 
compression of volume. The 4%, 6%, 8% and 10% corresponds 
to pressures: 1.7 GPa, 3.1 GPa, 4.7 GPa, 6.4 GPa, and 8.3 
GPa, respectively. 



However, it is necessary to point out that the shape of 
the Fermi surfaces is dependent on the magnetic moment 
which is sensitive to the exchange-correlation potential 
(LSDA or GGA) and the Fe-As bond length, and thus 
different calculations give quite different Fermi surfaces. 
Especially, the closed Fermi surfaces in the r — Y line, 
which are mainly d xy character, seem not to be a com- 
mon feature of the iron-based superconducting family. 
The cylinder tubes (see Fig. [5]), which are mainly d xz / yz 
character and show two-dimensional behavior under pres- 
sure, are the common features of the iron-based super- 
conducting family. Although the topological transition 
is observed under the pressure, such a transition is sensi- 
tively dependent on the amplitude of magnetic moment. 
As shown in Fig. [4l the two-dimensional cylinder tube 
in the T — Z line appears when the moment is about 1.3 
fiB in the GGA calculation. Considering that the exper- 
imentally observed magnetic moment is only about 0.8 
\xb-, the cylinder-tube Fermi surfaces of the striped an- 
tiferromagnetic state should appear at smaller pressure, 
because the current electronic structure calculations al- 
ways overestimate the magnetic moment but describe the 
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structural parameter very well. It should give some clues 
for the relationship between magnetism and supercon- 
ductivity, as discussed previously for LaFeAsO system.— 

Finally, wc investigate pressure effect on the band 
structures of the nonmagnetic state. As discussed above, 
the Fe-As bond length is always underestimated in the 
optimized structure of nonmagnetic state under pres- 
sure, consequently d xy bands are below Fermi level at T 
point (See Fig. [IJ, only shift downwards a little with in- 
creasing pressure. While d X z,yz bands shift downwards 
significantly under pressure, which across Fermi level 
as compression of volume over 8%, corresponding to a 
pressure value of 4.7 GPa in nonmagnetic calculations. 
Such a crossover will lead to a Fermi surface topologi- 
cal transition of the hole-type cylinder-tube closed at T 
point. However, it should be mentioned that the non- 
magnetic results might not be true. Because ARPES 
experiments always indicate the hole-type Fermi sur- 
faces around T connects with larger gap in the iron- 
pnictide materials family, which should be important to 
superconductivit y 2 ^ 26 . Thus the topological transition 
induced by pressure predicted by nonmagnetic calcula- 
tion would not coincide with the fact that T c reach max- 
imum nearby such a pressure. 

It has been observed that the superconducting states 
coexist with the striped magnetic stated, but at present 
we cannot say that the mechanism of the superconductiv- 
ity is based on the magnetism. It is just believed that the 
superconductivity in the Fe based materials is unlikely 
mediated by phonons, but likely originated from some 
magnetic or spin fluctuations. Our calculated results in- 
dicate that the GGA calculations yield reasonable lattice 
constants and za s only when we consider the magnetism. 
This implies that the magnetic effect is essential in these 
materials. Since the d xz / yz bands around Fermi level are 
determined by the Fe-As antibonding, as the distance of 
Fe-As increases, the k z related hoppings are weakened 
and the bands become narrow, the large experimental 
ZAs parameter leads to typical two-dimensional proper- 
ties. In current ab initio electron-phonon calculations, 
the optimized structure of nonmagnetic state is used, 
which results in the underestimation of the Fe-As bond 
length and more three-dimensional character of Fermi 
surfaces, that would influence properties of phonon and 
electron-phonon interaction. Such a magnetism related 
effect might have substantial influence on the electron- 
phonon coupling. A calculation including the magnetic 
effect are highly expected to investigate how much the 
magnetic background influences the electron-phonon in- 
teraction, to reveal the possibility of that phonon might 
play some important role on the superconductivity. On 
the other hand, LDA calculations yield smaller magnetic 
moment and underestimates the zas parameter. The cor- 
responding Fermi surfaces is substantially different from 
that of ARPES experiments^. If the superconductiv- 
ity is not directly dependent on the static moment, but 
determined by the lattice parameters and spin fluctua- 
tions, the changing of the moment value does not affect 



the superconductivity. 



V. COMPARED WITH EXPERIMENTAL 
RESULTS INVOLVED 

The main implication of pressure effect on the crys- 
tal structure of BaFe2As2 is to shorten the Fe-Fe bond 
length and the Fe-As bond length, but to maintain As- 
Fe-As bond angles simultaneously We could compare 
the pressure effect on crystal structure with the doping 
effect in (Bai_ a; K 2 ;)Fe2As2i. With the doping of K into 
BaFe2As2, T c goes up and reaches a maximum of 38 K 
at £=0.4. The impressive change is that As-Fe-As an- 
gle e becomes the ideal tetrahedral angle of 109.5° as 
x=0.4. Under pressure, the FeAs4 - tetrahedra is hard 
to be changed and still far away from the ideal structure. 
It seems to imply that superconductivity could not be 
improved simply by hydrostatic pressure in BaFe2As2. 

In the iron-pnictide superconducting family, the crystal 
structure, magnetism and superconductivity play com- 
plex game. It has been reported that a structural transi- 
tion of c-axis collapse occurs in CaFe2As 2 27l2S under pres- 
sure. Such a structural phase transition happens nearby 
the upper boundary of superconducting dome. In the 
striped antiferromagnetic calculations of BaFe2As2, we 
find a collapse of c-axis as the compression of volume 
above 8%, which should be relevant to the topological 
transition of Fermi surfaces. As the topological transition 
of Fermi surfaces influence the magnetic moment of iron 
slightly, the related c/a collapse should not be relative 
to a magnetic transition. In fact, we observe a magnetic 
collapse accompanying by a sudden softness of c/a ratio 
as lattice volume is compressed over 16%, corresponding 
to near 15 GPa pressure. It is far away from the known 
superconducting dome (2-6 GPa) up to date. We have 
performed LSDA calculations to examine the influence of 
overestimated magnetic moment by GGA. We found that 
the LSDA calculated magnetic moment is about 0.35 /is 
smaller than GGA, and the magnetic collapse happens at 
about 14% compression of volume, corresponding to pres- 
sure about 13 GPa. These results are consistent with re- 
cent experimental report that a suppression of magnetic 
order by pressure is observed as pressure above about 13 
GPa^. 

Now, we attempt to consider the implication of the 
striped antiferromagnetic and nonmagnetic calculations 
for paramagnetic state. It has been found experimen- 
tally that superconductivity coexist with static magnetic 
order under pressure for BaFe2As2^, so that a strong 
magnetic background is essential to superconductivity in 
these families. First of all, in order to describe chemi- 
cal bonding more accurately, the magnetic effect should 
be considered. Since ferromagnetism is not stable under 
a large range of pressure, only antiferromagnetic calcula- 
tions is possible. As the optimized internal parameter za s 
of the striped antiferromagnetic state is in good agree- 
ment with experiments, it might be reasonable to use 
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these parameters to describe the crystal structures and 
electronic structures in paramagnetic state than those 
nonmagnetic results. It coincides with the fact that no 
significant structural transformation is observed as sys- 
tem goes from antiferromagnetic state to superconduct- 
ing state or paramagnetic state, except the compression 
of b/a orthorhombic transition induced by the antiferro- 
magnetic interaction. If the optimized lattice parameters 
obtained from the striped antiferromagnetic calculations 
are adopted, it turn out that d xz , d yz and d xy bands 
always cross Fermi level except for shifting downwards 
a little under pressure. For such a treatment of para- 
magnetic state, we found that the energy splitting be- 
tween dxz.yz and d xy bands decreases from 0.06 eV to 
0.04 eV at T point as volume is compressed 6%, which 
indicates that the d xy bands are determined by direct 
Fe-Fe d xy interaction while the d xz / yz bands are mainly 
determined by Fe d X z,yz - As p x / y interactions. Neverthe- 
less, the band shapes are maintained under pressure, and 
the Fermi surfaces are quite similar to those calculations 
with experimental lattice parameters (see Fermi surfaces 
of left column in Fig. [T]), except the sheets around T 
point expand more. Thus the Fermi surfaces are robust 
against the pressure, no topological transition is found as 
the lattice volume compressed up to at least 10%. The 
impressive features is very different from those nonmag- 
netic calculations. 

It should be reminded that the band structures near 
Fermi level, as well as Fermi surfaces are sensitive to Fe- 
As bond length and As-Fe-As angle, which highly depend 
on the details of the spin-polarization. Our calculations 
indicate that the shape of bands and Fermi surfaces are 
almost decided by the value of large range magnetic mo- 
ment. Therefore, if the pressure induces a high spin po- 
larized state transform to a low spin state (even nonmag- 
netic state) , not only d xy bands shift down to cross Fermi 
level, but also d xz ^ yz bands cross Fermi level. This would 
certainly influence the coupling between different Fermi 
surfaces, then strongly influence the electronic proper- 
ties. It would be very interesting to examine the detail 
of magnetism under pressure by experiment furthermore. 



VI. CONCLUSION 

In summary, we have employed an accurate all- 
electronic full-potential linearized augmented plane wave 
method within the density-functional theory to investi- 
gate the effect of pressure on the crystal structure and 
electronic properties of BaFe2As2- The calculated lat- 
tice parameters of the striped antiferromagnetic state is 
in good agreement with experiment, while nonmagnetic 
state calculation underestimates c-axis and Fe-As bond 
lengths and yields a smaller za s - We find that the c-axis 
is compressed more easily than a-axis in nonmagnetic 
calculation, while the antiferromagnetic state calculation 
lead to a nearly isotropic compressibility with pressure 
up to 6 GPa. With increasing pressure, the FeAs4 tetra- 
hedra changes little. For the striped antiferromagnetic 
state, we find that the Fermi surfaces get closer and closer 
to the r point in the ±z directions when pressure in- 
creases from 2 GPa to 6 GPa, corresponding to the pres- 
sure region where superconductivity occurs. We observe 
a Fermi surface topological transition around T point as 
the pressure beyond 6 GPa. The transition should cor- 
respond to a small c-axis collapse. A magnetic collapse 
inducing the softening of c-axis is found when pressure 
reaches 13 GPa. These first-principles results should be 
useful to understanding the antifcrromagnctism and elec- 
tronic states in the FcAs-bascd materials, and may have 
some useful implications to the superconductivity. 
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